Quantum dynamics of a driven three-level Josephson-atom maser 
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Recently, a lasing effect has been observed in a superconducting nano-circuit where a Cooper pair 
box, acting as an artificial three-level atom, was coupled to a resonator^. Motivated by this exper- 
iment, we analyze the quantum dynamics of a three- level atom coupled to a quantum- mechanical 
resonator in the presence of a driving on the cavity within the framework of the Lindblad master 
equation. As a result, we have access to the dynamics of the atomic level populations and the photon 
number in the cavity as well as to the output spectrum. The results of our quantum approach agree 
with the experimental findings. The presence of a fluctuator in the circuit is also analyzed. Finally, 
we compare our results with those obtained within a semiclassical approximation. 
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I. INTRODUCTION 

In recent years, a considerable progress has been 
achieved in the field of quantum manipulation with nano- 
circuits based on the Josephson effect^^^. This progress 
has been initially inspired by the ideas of quantum in- 
formation processing^, which require the physical real- 
ization of qubits as an elementary quantum information 
unit. The degree of control achieved in Josephson-based 
qubits is so high that these systems have become a test- 
bed for the ideas of quantum mechanics such as quantum 
noise detection^^, quantum measurements^'^, and real- 
ization of circuit-quantum electrodynamics (QED)^^'^^. 
Very recently, the experimental realization of a single 
Josephson-atom laser has been reported^. 

In the experiment of Ref. 1, a Cooper pair box (CPB) is 
coupled to a superconducting waveguide resonator. The 
CPB is used as a three-level artificial atom. While the 
lowest two levels constitute a qubit, population inver- 
sion is achieved with the Josephson quasiparticle (JQP) 
cycle^^ involving the third level. The lasing condition 
can be determined from the steady state photon num- 
ber, pointing out that a too strong pumping suppresses 
the lasing action^^. Experimentally, evidence for lasing 
action was found through measurements of the output 
power spectrum of the resonator. An additional driv- 
ing was applied on the cavity to induce phase locking, 
thereby enhancing the lasing effect. The results are con- 
sistent with theoretical work on a two-level atom cou- 
pled to a resonator Studies concerning the coupling 
between a superconducting qutrit and a resonator have 
been also carried out^^. However, currently no quantita- 
tive results are available for the dynamics of a three-level 
Josephson atom coupled to a resonator. 

In this paper we present a theoretical analysis of 
the lasing effect observed in the experiment^, using a 
quantum-mechanical approach. We first obtain the com- 
plete time evolution of the system, including transient 
effects upon switching on the pumping. We obtain esti- 
mates for the characteristic time scales of the correspond- 



ing dynamics. Then, we study the output spectrum of 
the cavity field, with and without an additional driving 
applied to the cavity. We show that the latter requires 
a full quantum treatment and cannot be obtained from 
a semiclassical approximation. Our results are in good 
agreement with the experimental findings. The possi- 
bility to induce off-resonance lasing, observed in the ex- 
periment with a second hot spot, is also considered by 
adding a two level fluctuator in the circuit. Since our 
model is system-independent it can be applied to the 
study of other circuit-QED implementations, such as the 
currently much-studied transmoni^. 



II. MODEL 

The system under consideration is depicted in 
Fig. [TJ It is composed of the three-level artificial atom 
{|0), |2)}, with an energy difference huio between the 
ground state and the first excited state, and the cavity 
with a mode frequency ujo/27t. These two subsystems are 
coupled coherently according to the Jaynes-Cummings 
Hamiltonian 

H = ^huJiQCTz + hujoo^a + ihg{(jQia^ — aioa), (1) 

where (Jij — Oz = cru — aoo and a (a^) is the canon- 

ical bosonic annihilation (creation) operator of a photon 
in the cavity. The dynamics of the third level |2) is de- 
scribed with a Lindbladian, as presented below. This 
Hamiltonian is obtained after applying the rotating wave 
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FIG. 1. Scheme of the three- level system coupled to a cavity 
with the corresponding transition rates. 
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FIG. 2. (Color online) Time evolution of the photon number 
in the cavity. The coupling and various rates correspond to 
the experimental parameters at the resonance, as discussed 
in the text. Inset: Distribution of the photon population in 
the steady state (histogram) compared to the corresponding 
binomial distribution (solid line). 

approximation (RWA), valid when the coupling strength 
is small compared to the typical frequency of the isolated 
subsystems, which will be the case in the following. The 
legitimacy of the RWA is also established through nu- 
merical simulations. The atom is pumped from |0) to |2) 
at the rate F, the state |2) decays to |1) at the rate 721. 
The reverse processes occur at the rates, respectively, 
720 and 712. Finally, the cavity has a damping rate a^, 
conferring to the photons a lifetime K:~^. As a conse- 
quence, the time evolution of the total density matrix p 
has two contributions: the evolution due to the coherent 
coupling between the artificial atom and the cavity ac- 
cording to Hamiltonian ^ and the evolution due to the 
incoherent processes and controlled by the Lindbladian 
L = Lr ^ I/721 + ^720 + ^712 + ^/^ with 

for the pumping and the relaxation rates, noting F = 702, 
and 

L^p = ^K{2apa^ — a^ap — pa^a)^ (3) 

for the damping of the cavity. These expressions are ob- 
tained in the Born-Markov approximation, supposing a 
weak coupling between the system and the environment. 
The resulting time evolution of the density matrix satis- 
fies the master equation^^- 

m = ^[H,pit)]+Lpit). (4) 

To characterize the coherence properties of the emitted 
field we calculate the output spectrum S{uj) defined as 
the Fourier transform of the cavity correlator 

/+00 
dre-^^^(a^(t + r)a(t)). (5) 




FIG. 3. (Color online) Time evolution of the level populations 
at the resonance (level |0) in blue, |1) in green, and |2) in red). 
The dynamics of the first pump cycles for a qubit initially in 
the state |0) is presented in the inset. 

The output spectrum is obtained from the steady-state 
density matrix using the quantum regression theorem^^, 
valid within our approach, which establishes a matrix 
relation between them. We also include the possibility 
to drive the cavity with the additional pumping 

Hd = ihv{e-''^^o) - e^^^a), (6) 

where w is the detuning from the cavity frequency. The 
amplitude can be expressed in terms of the photon num- 
ber A^o created by the driving^^ v = If the 
emitted field is coherent, the injection locking effect oc- 
curs and the cavity field oscillates at the same frequency 
as the driving field uoq ^ w. 

III. CHOICE OF PARAMETERS 

The results that we present below have been obtained 
using a particular choice for the numerical values of the 
various system parameters, corresponding to those of the 
experiment^. The three- level atom is a CPB, the proper- 
ties of which are controlled by an external parameter, the 
dimensionless gate- voltage Ug. Varying the parameter Ug 
corresponds to rotating the charge basis around the state 
|2) by an angle defined by tan2l9 = Ej/[Ec{ng - 1)], 
where Ej and Ec are the Josephson energy and the 
charging energy of the CPB, respectively {Ec/Ej 15). 
The qubit energy then reads huio = Ej/ sin26>, the cou- 
pling varies like sin26>, the rates F and 721 are propor- 
tional to cos^ ^, and the rates 720 and 712 are propor- 
tional to sin^ 0. In the experiment, the cavity frequency 
is cJo/27r 10 GHz; the resonance condition cjio = ujq 
for the lowest two atom levels and the cavity is achieved 
when the parameter = 1.1. At this working point, 
the atom-cavity coupling frequency is g/27T = 44 MHz. 
Population inversion is achieved using the JQP cycle; 
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FIG. 4. (Color online) Density plot of the spectrum (loga- 
rithmic scale) from our fully-quantum model as a function of 
the probing frequency and the reduced gate charge. 



the relevant rates are F = 4.2 GHz, 721 = 3.3 GHz, 
720 = 0.29 GHz, and 712 = 0.37 GHz. The damping rate 
is set to = 8.2 MHz. What can be measured experi- 
mentally is the spectrum of the cavity, Eq. (|5]), with or 
without an additional driving Eq. (|6j). The photons be- 
ing emitted at an energy of 10 GHz, this lasing effect is 
actually a masing effect. 



IV. DISCUSSION OF THE RESULTS 

The Lindblad master equation Eq. dU gives access to 
the time evolution of the photon number in the cavity 
and of the Josephson atom level populations. The dy- 
namics of these quantities is shown in Figs. [2] and [H using 
the parameters given above. We express time in units of 
the inverse pumping rate 1/F. We see that the transient 
time, i.e., the time needed to reach the steady state, is 
on the order of 4000 pumping cycles (a microsecond for 
the experiment). At very short time scales, the three- 
level atom shows a significant population imbalance; this 
accompanies a fast increase in the photon number in the 
cavity. In the steady state more than one hundred pho- 
tons are present in the cavity, in agreement with the ex- 
perimental estimates. The photon distribution follows a 
binomial law (see Fig. O inset), characteristic for corre- 
lated particles at zero temperature^. We next calculate 
the output spectrum Eq. ([5j), as a function of the fre- 
quency uj and of the parameter Ug. The result is plotted 
in Fig. m It presents a peak centered at the resonance. 
Furthermore, in the experiment the presence of charge 
fluctuations widens the spectrum. This broadening can 
be overcome by driving the cavity. Figure [5] represents 
the spectrum as a function of the driving strength. As v 
increases, the initial Lorentzian is converted into a Dirac 
peak located at the driving frequency, thus emphasizing 
the lasing effect. Charge noise can be strongly suppressed 
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FIG. 5. (Color online) Spectrum in the presence of an addi- 
tional driving on the cavity {vj/2ti — —1 MHz). The damping 
n has been increased fivefold. 



if a transmon qutrit is used instead of the standard CPB. 



V. COUPLING TO A TWO-LEVEL SYSTEM 

The experimental spectrum in the (ug^uj) plane shows 
an additional peak centered at ujiq ^ 1.5 ujq. Such a 
second peak is absent in our simulations, which is con- 
sistent with the fact that the coupling strength g is be- 
low the threshold for two-photon masing. Indeed, even 
when including the counter-rotating terms of the master 
equation that are neglected in the RWA and can lead 
to multi-photon processes, lasing effect occurs only at 
the resonance. Moreover, if the observed peak were due 
to two-photon processes, it should have been located at 
ujiQ c^2ujo. A possible source of additional resonances is 
the coupling to a two- level system (TLS). We will con- 
sider two kinds of coupling, first a resonant coupling with 
the lasing transition of the CPB and second a dispersive 
coupling with both the CPB and the cavity. To under- 
stand if it is possible to observe off-resonant lasing effect 
with a simple TLS we focus on the steady state photon 
number. 

The Josephson junction of the CPB can be a source 
of fluctuators, resulting from the tunneling of a charge 
between two sites in the insulating layer^^. The TLS is 
then composed of two states, ground state \g) and excited 

state |e), separated by an energy hwf = y£^+4T^, 

where Ef is the energy difference between the two sites 
and T is the tunneling strength. The TLS is described by 
the Hamiltonian Hf = ^huufc^z, with = [^+7^-] where 
= \e){g\ and ^_ = The tunneling charge position 
couples to the Cooper pair number, with an energy %r, 
forming a four level system {|0), |1)} {|^), |e)}. This 
system is furthermore coupled to the single electron state 
1 2) with the incoherent pumping and to the cavity states 
through the Jaynes-Cummings Hamiltonian ([1]). While 
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FIG. 6. (Color online) Steady state photon number in the 
presence of a two-level system. The different lines correspond 
to different coupling strengths in the case of a resonant cou- 
pling and a dispersive one (see the legend, g/2n = 44 MHz). 
The inset is a zoom around the frequency 1.5 cjq, where the 
second peak appears. The damping n has been increased five- 
fold and the photon number is normalized by the number at 
the resonance without TLS. 
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FIG. 7. (Color online) Maximum value of the spectrum 
as a function of the reduced gate charge. The quantum 
solution is plotted in purple and the semi-classical one in 
black. Inset: Accuracy of the factorization Eq. (|10|) . The 
time evolution of the real part of the normalized difference 
[{sz)S{t) — C{t)]/ {(TzTi) in the rotating frame of the cavity is 
plotted at the resonance (cl;io = cjq) in blue and at the second 
peak of 5'max(Tip) (cJio = 1.06 cjq) in red. 



varying the gate voltage, different transitions of the four 
level system can become resonant with the cavity and 
induce lasing. The interaction Hamiltonian turns out to 
comprise both a transverse and a longitudinal coupling, 
as well as frequency shifts, 

(7) 

where viq = 4.grEc{l - ng)/huJio, Uf = QrEfil - 
ng)/hujf, gl = -grEjT/h^ujiQUjf, and g[ = 
2gr{l — ng)EcEf /Ji^ujioujf. The energy spectrum of 
the four level system {iV^i), |'02), IV^a), |'04)} is given 
by £^1,4 tK^io + + %z and £^2,3 = 

T^y^(^io^^^/)M^4^/2 - hgi, where wioj = ujioj + 
z^io,/. The ground state is = |0,^), the high- 

est state is |V^4) = |l,e) and the central terms are ob- 
tained after rotating |0, e), |1,^) by an angle a satisfying 
tan 2a = 2gt / {vo f — tuio). Finally, the coupling Hamilto- 
nian with the cavity, obtained from Eq. ([T]), reads 

^FLS = % [cos a(6'o2 + ^13) + sin ^(^oi - ^'23)] + H.c, 

(8) 

where Sij = A lasing effect thus occurs if the 

transition 2-0, 3-1, 1-0, or 3-2 is in resonance with the 
cavity and the corresponding coupling strength is large 
enough. The photon number as a function of the fre- 
quency ujio is plotted in Fig. [6] for different values of 
the coupling strength gr. The frequency of the TLS is 
adjusted close to ujq (±10%) to observe the second reso- 
nance at ujio = 1.5 cJo, in the regime Ef = T. 

A second resonance can also be produced by a fluctua- 
tor of frequency ujf = 1.5 cjq coupled to both the Joseph- 
son atom and the cavity. The cavity can then be in- 
directly excited when the transition c^io approaches the 



resonance with the TLS, uio — ojf. In this dispersive 
regime, we consider only the transverse coupling^—, 

Hd = h [g\{(J- + a+) + g\{a + a^)] (c^+ + c^_), (9) 

where g\i2 is the transverse coupling strength between 
the TLS and the Josephson atom/cavity. The coupling 
g\ has the same dependence on Ug as g while g^ is Ug- 
independent. The steady state photon number is plotted 
in Fig. [6] as a function of the gate voltage through uoiq for 
gl = g and different values of ^2/^- The damping rate 
has been increased fivefold for numerical reasons, but this 
does not change the results qualitatively. 

In both cases the effect of the TLS on the photon num- 
ber is small even for strong couplings (see Fig. [6l inset). 
An ultra-strong coupling between the fiuctuator and the 
system, unrealistic for the experimentii, is needed to ob- 
serve off-resonant lasing. A simple TLS is thus unlikely 
to explain the second hot spot of the experiment. 

VI. SEMI-CLASSICAL APPROXIMATION 

The results presented so far were obtained by numeri- 
cally solving the Lindblad master equation. When the 
steady-state photon number is large, one can use the 
semi-classical approximation to get analytical results 
It consists of factorizing the operators pertaining to the 
three-level atom and to the cavity. The steady-state 
value of the photon number, for instance, is then obtained 
from a set of equations involving (cr^), (coi^^^), and 

(cr^n). The latter can be factorized in the semi-classical 
limit and the resulting solution is in good agreement with 
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the numerical results. At the level of the spectral func- 
tion, the time derivative of S{r) induces more complex 
correlators such as C{r) = {az{r)a\r)a{0)). Using the 
amplitude-phase representation of the operator a and as- 
suming that the correlation time of the phase fluctuations 
is much longer than that of the amplitude fluctuations, 
the factorization can be improved^^ 

(a,(r)at(r)a(0)) ^ ^ (^(a,) + ^) 5(r). (10) 

The set of differential equations leads to a Lorentzian 
spectrum of width k = k/2 - 2g'^{T + 72i)(s^)/[(r + 
721)^ + 4(5^], and centered at the reduced frequency 
5uj = (4^^J(5;,))/[(r + 72i)^ +4(5^], where we note (s^) = 
\ {{^z) + 5 = uio — uq is the detuning. 

The maximum spectrum S'max with respect to the fre- 
quency uj is plotted as a function of the reduced charge 
gate Ug in Fig. [71 and compared to the Lorentzian solu- 
tion in the semi-classical limit. The time evolution of the 
correlators in Eq. (fTQ|) in the rotating frame of the cavity 
is shown in the inset, where the real part of the normal- 
ized difference [{sz)S{r) — C{r)]/ (a^n) is plotted for two 
different values of the gate voltage. At the resonance, the 
dynamics of the factorized correlator {sz)S{r) is in good 
agreement with C{r). Off-resonance at ujio = 1.06 cjo, 
where the semi-classical spectrum exhibits the second 
peak, the difference oscillates at ~ a;o/27r + l MHz, giving 
rise to a non-negligible contribution in the Fourier trans- 
form. These comparisons reveal that the semi-classical 
treatment is not correct in the region close to the res- 
onance where the correlations between the atom and 
the cavity cannot be neglected. The resulting double 



peak structure thus appears to be an artifact of the fac- 
torization^^. Further improvement of the factorization, 
Eq. (fTQ|) . is needed to describe the spectrum properly in 
the semi-classical limit. 



VII. CONCLUSION 

In conclusion, the Lindblad master equation together 
with the quantum regression theorem are powerful tools 
to calculate quantum mechanically the time-evolution of 
the photon number and the output spectrum of the cav- 
ity. A comparison with the experimental results of Ref. [l| 
gives access to the typical time-scales of the system. The 
calculation of the spectrum enables us to understand the 
experimental results and the effect of driving the cavity. 
It shows in particular that the presence of charge noise 
reduces the lasing effect. Considering the presence of a 
fluctuator in the system, we show that an ultra-strong 
coupling is needed to explain the second hot spot. Fi- 
nally, the fully quantum treatment for a three-level arti- 
ficial atom, based on the density matrix of the whole sys- 
tem, allows to figure out the validity of the semi-classical 
approximations, which do not take into account all the 
correlations between the atom and the cavity. 
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